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L  Introduction 


In  systems  where  independent  users  transmit  through  a  single  common  channel,  the  deploy¬ 
ment  of  random  access  transmission  algorithms  is  frequently  desirable,  for  the  following  rea¬ 
sons:  (1)  They  are  implemented  independently  by  each  user,  without  a  priori  coordination 
among  the  users.  (2)  They  are  insensitive  to  changing  user  population.  (3)  They  induce  low 
delays  for  low  input  rates,  when  the  user  traffic  is  bursty. 

In  this  paper,  we  propose  and  analyze  a  full  feedback  sensing  window  random  access  algo¬ 
rithm.  The  algorithm  was  first  proposed  for  systems  with  strict  delay  limitations,  [31,  and  it 
requires  that  each  user  know  the  overall  feedback  history,  (full  feedback  sensing).  As  compared 
to  other  such  existing  algorithms,  the  proposed  algorithm  has  the  following  interesting  proper¬ 
ties:  (1)  It  can  be  easily  modified  to  operate  under  limited  feedback  sensing,  where  each  user 
follows  the  feedback  history  from  the  time  he  generates  a  message  to  the  time  when  this  message 
is  successfully  transmitted.  (2)  When  the  Poisson  user  model  is  adopted,  the  algorithm  attains  the 
same  throughput  as  that  attained  by  the  Capetanakis’s  dynamic  algorithm,  [1],  while  it  induces 
lower  delays  for  arrival  rates  above  0.30,  and  better  resistance  to  feedback  channel  errors.  (3) 
The  simple  operations  of  the  algorithm  allow  analytical  evaluation  when  strict  delay  limitations 
exist,  [3].  Its  simplicity,  together  with  its  regenerative  properties,  provide  the  means  for  the 
analytical  evaluation  of  the  output  traffic  interdeparture  distribution  induced  by  the  algorithm. 
The  analysis  of  the  latter  distribution  is  important  when  several  systems  which  use  some  Ran¬ 
dom  Access  Algorithm,  (RAA),  for  internal  transmissions  interact,  and  it  is  not  quite  feasible 
when  either  the  Capetanakis’s,  [1],  or  the  Gallager’s,  [2],  algorithms  are  deployed.  (4)  As  com¬ 
pared  to  Gallager’s  algorithm,  [2],  the  proposed  algorithm  operates  in  systems  where  the  Poisson 
user  model  is  not  valid,  (e.g.,  when  more  than  one  packets  can  be  generated  within  a  given  time 
instant),  and  can  be  then  analyzed. 

The  organization  of  the  paper  is  as  follows:  In  section  II,  the  system  model  is  presented, 
and  the  algorithm  is  described.  In  section  III,  throughput  and  delay  analysis  for  the  Poisson  user 
model  is  included.  In  section  IV,  the  performance  of  the  algorithm  in  the  presence  of  feedback 
channel  errors  and  its  operation  under  limited  feedback  sensing  are  discussed.  In  section  V,  the 
output  traffic  interdeparture  distribution  induced  by  the  algorithm  is  analytically  evaluated.  In 
section  VI,  some  conclusions  are  drawn. 

n.  The  System  Model  and  the  Algorithm 

We  assume  packet  transmitting  users,  slotted  channel,  binary  collision  versus  noncollision, 
(CNC),  feedback  per  slot,  and  zero  propagation  delays,  and  absence  of  feedback  errors.  We  also 
assume  that  collided  packets  are  fully  destroyed  and  retransmission  is  then  necessary.  Time  is 
measured  in  slot  units,  slot  t  occupies  the  time  interval  [t,  t+1),  and  x,  denotes  the  feedback  that 
corresponds  to  slot  t;  x,  =C  and  x,  =  NC  represent  then  collision  and  noncollision  slot  t,  respec¬ 
tively.  For  this  system,  let  the  following  full  feedback  sensing  synchronous  random  access  algo¬ 
rithm  be  deployed. 

The  algorithm  utilizes  a  window  of  length  d.  Let  t  be  a  time  instant  that  corresponds  to  the 
beginning  of  a  slot  such  that,  for  some  t,  <  r,  all  the  packet  arrivals  in  (0,  t  j]  have  been  success¬ 
fully  transmitted  by  the  algorithm  and  there  is  no  information  regarding  the  arrival  interval 
(r i ,  r).  The  instant  t  is  then  called  Collision  Resolution  Point,  (CRP),  the  arrival  interval  (0,  f|]  is 
called  "resolved  interval",  and  the  interval  fri,  t]  is  called  "the  lag  at  t".  In  slot  t,  the  packet 


arrivals  in  (rl4  /j^nrinf/!  +  A,fj]  attempt  transmission,  and  the  arrival  interval  (t lt  t2]  is  then 
called  the  "examined  interval".  If  (t\,  t2\  contains  at  most  one  packet,  then  it  is  resolved  at  t.  If 
fr  i ,  t2]  contains  at  least  two  packets,  instead,  then  x,  =  C,  a  collision  occurs  at  t,  and  its  resolution 
starts  with  slot  t+1.  Until  the  collision  at  t  is  resolved,  no  packets  that  have  arrived  after  t2  are 
allowed  transmission.  The  time  period  required  for  the  resolution  of  the  latter  collision  is  called 
the  Collision  Resolution  Interval,  (CRI).  If  the  examined  interval  contains  at  most  one  packet, 
then  the  CRI  lasts  one  slot.  During  some  CRI  which  starts  with  a  collision  slot,  each  user  acts 
independently  via  the  utilization  of  a  counter  whose  value  at  time  t  is  denoted  rt.  When  a  user 
transmits  for  the  first  time  he  sets  r,  =  1  .  The  counter  values  can  be  either  1  or  2,  and  they  are 
updated  and  utilized  according  to  the  rules  below. 

1.  The  user  transmits  in  slot  t,  if  and  only  if  r,  =  1.  A  packet  is  successfully  transmitted  in  t,  if 
and  only  if  r,  =  1  and  x,  =  NC. 

2.  The  counter  values  transition  in  time  as  follows: 

(a)  If  x,_i  ss  NC  and  =  2,  then  r,  =  1 

(b)  If  x,_i  =  C  and  r,^  =  2,  then  r,=2 

(c)  If  x,-i  *  C  and  r,_ j  =  1,  then 

J 1,  w.p.  0.5 
r'  =|  2,  w.p.  0.5 

From  the  above  rules  it  can  be  seen  that  a  CRI  which  starts  with  a  collision  slot  ends  with  two 
consecutive  noncollision  (NC)  slots.  Furthermore,  two  consecutive  NC  slots  can  not  occur  at 
any  other  instant  during  a  CRI.  Thus,  the  end  of  a  such  a  CRI  can  be  identified  by  all  the  users  in 
the  system,  upon  the  observation  of  two  consecutive  NC  slots. 


Remarks  We  note  that  the  algorithmic  operations  can  be  depicted  by  a  two-cell  stack,  where  at 
each  time  instant  t,  cell  1  contains  the  transmitting  users,  (those  with  r,  =  1),  and  cell  2  contains 
the  withholding  users,  (those  with  r,  =  2).  The  algorithm  lumps,  thus,  the  unsuccessful  users 
together.  In  contrast,  Capetanakis’s  algorithm  distributes  the  unsuccessful  users  across  the  cells 
of  an  infinite-cell  stack.  As  with  Capetanakis’s  dynamic  algorithm,  the  window  size  A  is  here 
subject  to  optimization  for  throughput  maximization. 

III.  Throughput  and  Delay  Analysis 

In  this  section,  we  present  throughput  and  delay  analysis  of  the  algorithm.  We  assume  con¬ 
tinuous  feedback  sensing  and  we  adopt  the  Poisson  user  model.  Indeed  as  proven  in  [5],  for  a 
large  class  of  random  access  algorithms  (RAAs),  as  the  population  size  increases  the  stability  of 
an  algorithm  in  the  class  is  determined  by  its  throughput  under  the  Poisson  user  model. 

Consider  the  system  model  and  the  algorithm  in  Section  II.  Let  the  system  start  operating  at 
time  zero.  Let  r, ;  /  >  1  be  the  sequence  of  successive  CRPs,  and  let  X,  be  the  lag  at  r*.  The 


sequence  X ;  i  £  1  is  a  Markov  Chain  with  state  space  F,  F  is  an  at  most  denumerable  subset  of 
the  interval  [1,  <» ).  It  can  be  seen  that  any  state  can  be  reached  from  any  other;  therefore  X, ;  /'  >  1 
is  an  irreducible  Markov  Chain.  Since  P  (X1+]  =  1 1  X,  =  1)  >  0  ,  we  conclude  that  Xt\  i  >  l  is  also 
aperiodic.  Thus  Pake’s  Lemma  [11]  applies,  and  gives  that  the  following  condition  is  sufficient 
for  the  ergodicity  of  the  Markov  Chain,  (stability  of  the  system): 

£(/  I  A,d)  <  A  (1) 

where  £(/  I  A ,d)  denotes  the  expected  length  of  a  CRI,  given  that  it  starts  with  an  examined 
interval  of  length  A  and  with  a  lag  of  length  d. 

Since  the  Markov  Chain  is  uniformly  downward  bounded  (i.e.  there  exists  a  constant  m 
such  that  the  transition  probabilities  pkj  satisfy  pkj- 0  for  j  <k-m.  Here  m=  A  -  1),  Kaplan’s 
Theorem  [12]  applies  and  gives  that: 

If 


£(/  IA,d)  >  A 

then  the  Markov  Chain  is  not  ergodic  and  the  system  is  unstable. 

Let  Lk  denote  the  expected  length  of  a  CRI  given  that  it  starts  with  a  collision  of  multipli¬ 
city  k.  We  can  then  write: 


£(/  IA  ,d)=  ££(/  IA  M  =  ZLke-**^£- 


** o 


it! 


*= o 


it! 


since 


(2) 


(2a) 


£(/  IA,rf,it)  =  L* 

depends  only  on  k. 

In  the  Appendix  we  show  that: 

(i)  Lk ;  k  >  0  can  be  computed  recursively,  and 

(ii)  Lk  are  quadratically  upper  bounded,  Lk<Li  =  ak2  +$k+y,  k  2  2. 

Expression  (1)  together  with  (i)  and  (ii)  are  used  in  the  computation  of  the  algorithmic 
throughput 

We  define  the  delay  £>„,  experienced  by  the  n-th  packet  as  the  time  difference  between  its 
arrival  and  the  end  of  its  successful  transmission.  We  are  interested  in  evaluating  the  first 
moment  of  the  steady  state  delay  process,  when  it  exists.  Let  r,  =  1  ,  X\  =  1  ,  and  define  TM  as 
the  first  CRP  after  T,  at  which  the  lag  has  length  one.  From  the  description  of  the  algorithm  it 
can  be  seen  that  the  induced  delay  process  probabilistically  restarts  itself  at  the  beginning  of 
each  slot  T, ,  i  =  1 ,2  , ...  The  interval  [7'1,7’I+I)  will  be  referred  to  as  the  i-th  session.  Note  that  the 
sessions  have  lengths  that  are  i.i.d  random  variables. 

Let  /?,;/  =  1  ,2 ... ,  denote  the  number  of  packets  successfully  transmitted  in  the  interval  (0, 
£,  ]  ;  (note  that  /?,  also  represents  the  number  of  packets  arrived  during  the  interval  [0,  T,  -  1), 
since  r,  is  a  CRP  at  which  the  lag  has  length  one).  Then,  Q,  =  /?1+1  -/?, ;  i  >  1,  is  the  number  of 
packets  successfully  transmitted  in  the  interval  ( r„rl+i  ] ,  these  are  the  packets  that  arrived  dur¬ 
ing  the  interval  [  T,  -  l,r(+1  -  l ).  The  sequence  R, ;  i  £  1,  is  a  renewal  process  since  (2, ;  i  1 1,  is  a 
sequence  of  nonnegative  i.i.d  random  variables.  Furthermore,  the  delay  process  Dn  ,n>  1,  is 

regenerative  with  respect  to  the  renewal  process  /?,  ;/£  1 ,  with  regeneration  cycle  Q\.  From  the 

Q\ 

regenerative  theorem  [4],  we  conclude  that  if  Q  =  E(Q\)  <  «  and  W  =  £  {  £  DJ  <  <*>,  then  there 


(3) 


exists  a  real  number  D  such  that, 

D  =  lira /i_1  =  Iim/J-I£(£D,)  =  —  w.p.  1 

i-i  2 

In  addition,  since  P  (Q\  =  1)  >  0,  the  distribution  of  Q  j  is  aperiodic  and  there  exists  a  random 
variable  such  that  the  sequence  D*  ;  n  >  1  converges  in  distribution  to  D„  represents  the 
steady  state  delay  induced  by  the  algorithm  and  its  mean  satisfies  the  equality 

E(D„)  =  ~  (4) 

The  quantity  D  will  be  referred  to  as  the  mean  packet  delay.  From  (4)  we  observe  that  the  mean 
packet  delay  can  be  determined  by  computing  the  quantities  of  the  right  hand  side  of  the  equal¬ 
ity.  In  the  Appendix  we  develop  two  systems  of  linear  equations  whose  solution  may  be  used  to 
compute  the  mean  cycle  length  Q  and  the  mean  cumulative  delay  W. 

In  Table  1,  we  include  the  computed  upper  and  lower  bounds,  D“  and  D*  respectively,  on 
the  mean  packet  delay  D,  for  various  Poisson  intensities  X,  and  for  both  the  proposed  and  the 
Capetanakis’s  dynamic  algorithms. 


Proposed 

algorithm 

Capetanakis  dynamic 
algorithm 

X 

Dl 

Du 

Dl 

Du 

ipm 

1.562 

1.563 

1.563 

1.564 

EES 

1.708 

1.716 

1.713 

1.719 

0.10 

1.888 

1.917 

1.903 

1.921 

0.16 

2.257 

2.363 

2.308 

2.362 

2.607 

2.812 

2.712 

ItPI 

3.103 

3.467 

3.308 

3.476 

0.30 

4.412 

5.197 

4.976 

5.365 

0.32 

5.162 

6.170 

5.973 

6.501 

0.36 

7.941 

9.665 

9.798 

10.883 

0.38 

11.008 

14.121 

0.40 

18.262 

24.427 

0.42 

57.354 

67.665 

78.530 

90.212 

Table  1 


Upper  and  Lower  Bounds  on  Steady-State  Expected  Delays 


Regarding  the  throughput  X*  and  the  optimal  window  size  A*,  the  following  results  were  found. 


Proposed  Algorithm: 

X’  =  0.4295 

A*  =  2.33 

Capetanakis’  Dynamic  Algorithm: 

Table  2 


Throughputs  and  Optimal  Window  Sizes 

From  Table  2,  we  observe  that,  for  the  Poisson  user  model,  the  algorithm  in  this  paper 
attains  the  same  throughput  as  the  Capetanakis’s  dynamic  algorithm,  but  utilizes  a  smaller  win¬ 
dow  size.  From  Table  1,  we  observe  that  the  two  algorithms  induce  practically  identical  delays 
for  Poisson  rates  in  (0,  0.30),  while  for  Poisson  rates  in  (0.30,  0.42],  the  proposed  algorithm 
induces  lower  delays. 

Remarks  It  may  seem  surprising  that  the  algorithm  in  this  paper  attains  the  same  throughput  as 
the  Capetanakis’s  dynamic  algorithm,  and  that  it  outperforms  the  latter  in  terms  of  delay  perfor¬ 
mance.  Indeed,  the  expected  lengths  Lk  in  (2a)  are  bounded  by  quadratic  expressions,  while  the 
same  lengths,  Lk,  for  the  Capetanakis  algorithm  are  bounded  by  linear  functions  of  k.  However, 
Lj  =4.5  while  L\  =  5.  At  the  same  time,  since  A*  =  2.33  for  the  proposed  algorithm,  the  probabil¬ 
ity  of  a  higher  than  two  multiplicity  collision  is  very  small.  The  multiplicity-two  events  thus 
prevail,  and  the  algorithm  in  this  paper  becomes  better  than  the  Capetanakis’s  algorithm.  We 
note  that,  as  found  in  [3],  the  algorithm  performs  very  well  in  environments  where  strict  delay 
limitations  exist.  Then,  it  allows  significant  improvement  in  delay  performance,  at  the  expense 
of  minimal  traffic  loss.  In  addition,  the  analysis  of  the  algorithmic  performance  when  strict 
delay  limitations  exist  is  relatively  simple,  while  the  same  analysis  for  the  algorithms  in  [1]  and 
[2]  is  still  an  open  problem. 


IV.  Performance  in  the  Presence  of  Feedback  Errors  and  Operations  in  the  Limited  Sens¬ 
ing  Environment 

In  this  section,  we  study  two  important  characteristics  of  the  algorithm.  Namely,  its  perfor¬ 
mance  in  the  presence  of  feedback  channel  errors,  and  its  operation  and  performance  under  lim¬ 
ited  feedback  sensing. 

Performance  in  the  Presence  of  Feedback  Errors 

Let  us  assume  that  due  to  noisy  conditions,  the  following  types  of  feedback  errors  may 
occur:  With  probability  e  an  empty  slot  may  be  seen  by  the  users  as  a  collision  slot.  Also,  with 
probability  8  a  slot  occupied  by  a  single  transmission  may  be  seen  by  the  users  as  a  collision  slot. 
Let  us  also  assume  that  a  collision  slot  is  always  recognized  correctly  by  the  users.  We  consider 
the  case  where  the  probabilities  e  and  5,  being  system  characteristics,  are  known  a  priori.  Then 
given  e  and  5,  the  window  size  A  is  optimized  for  throughput  maximization.  We  performed 
throughput  analysis,  (the  details  are  included  in  the  Appendix),  for  both  the  proposed  algorithm 


and  the  Capetanakis’s  dynamic  algorithm,  [1].  We  exhibit  the  results  in  Table  3.  From  Table  3, 
we  conclude  that  the  proposed  algorithm  is  very  insensitive  to  feedback  errors.  Even  for  the 
practically  extreme  case  e  =  5  =  0.1,  the  throughput  is  almost  90%  of  its  value  in  the  error  free 
case.  We  also  conclude  that  the  proposed  algorithm  allows  operation  (positive  throughput)  as 
long  as  e  <  1  and  5  <  1,  while  if  eS0.5  the  throughput  for  the  Capetanakis’s  algorithm  is  then 
zero.  We  notice  for  example  that  for  e  =  0.5  and  5  =  0,  the  proposed  algorithm  attains  throughput 
as  high  as  X'  =0.325. 


e 

5 

X  proposed  alg. 

X '  Cap. 

0.00 

0.00 

0.4295 

0.00 

0.01 

0.4248 

0.4258 

0.00 

0.3873 

0.00 

0.3463 

:'U 

0.00 

0.2655 

0.2731 

0.00 

0.50 

0.2251 

0.2310 

0.01 

0.00 

0.4272 

0.4272 

0.10 

0.00 

0.4117 

0.4043 

0.20 

0.00 

0.3930 

0.40 

0.00 

0.3503 

0.45 

0.00 

0.3382 

0.50 

0.00 

0.3503 

0.0000 

0.60 

0.00 

0.3050 

0.0000 

0.70 

0.00 

0.2750 

0.0000 

0.80 

0.00 

0.2280 

0.0000 

0.90 

0.00 

0.1700 

0.0000 

0.10 

0.10 

0.3706 

0.3672 

0.20 

0.20 

0.3139 

0.2972 

0.30 

0.30 

0.2589 

0.2166 

0.40 

0.40 

0.2064 

0.1205 

0.30 

0.50 

0.1885 

0.1511 

0.30 

0.70 

0.1183 

0.0886 

0.90 

0.90 

0.0105 

0.0000 

Table  3 

Throughputs  as  a  function  of  e  and  5.  Window  sizes  optimized  for  every  pair  (e  ,  5) 


In  some  systems,  the  probabilities  e  and  6  may  not  be  known  a  priori.  In  this  case,  an  algorithm 
may  be  designed  subject  to  the  assumption  of  error  free  feedback.  We  found  the  throughputs  for 
the  proposed  and  the  Capetanakis’s  dynamic  algorithms,  in  this  case,  for  various  values  of  the 
error  probabilities;  that  is,  the  maximum  Poisson  rates  for  which  (1)  is  satisfied  subject  to  the 
constraint  that  the  windows  are  those  in  Table  2.  Our  results  are  shown  in  Table  4. 


X  prop.  alg. 


X  Cap. 


0.00 

0.00 

0.00 

0.00 

0.01 

0.10 

0.20 

0.3 

0.3305 

0.37 

0.3787 

0.10 

0.10 

0.10 

0.20 

0.20 

0.20 

0.30 

0.30 

0.30 


0.427 

0.408 

0.377 

0.320 

0.289 

0.214 

0.000 

0.363 

0.318 

0.273 

0.328 

0.279 

0.232 

0.261 

0.205 

0.153 


.429 

.410 

.391 

.262 

0.427 

0.401 

0.355 

0.248 

0.000 

0.000 

0.000 

0.361 

0.318 

0.272 

0.311 

0.266 

0.221 

0.192 

0.138 

0.088 


Table  4 


Throughputs  as  a  function  of  e  and  5.  Window  sizes  unchanged  for  every  pair  (e  ,  6) 


From  Table  4,  we  observe  that  the  proposed  algorithm  is  better  than  the  Capetanakis’s  algorithm. 
For  example  the  maximum  e  value  for  which  the  proposed  algorithm  is  stable  is  0.378,  while  the 
Capetanakis’s  algorithm  becomes  unstable  for  e  £  0.33. 


Operations  in  the  Limited  Sensing  Environment 


In  the  limited  sensing  environment,  it  is  required  that  each  user  monitor  the  channel  feed¬ 
back  only  from  the  time  he  generates  a  packet,  to  the  time  this  packet  is  successfully  transmitted. 
Therefore,  the  users’  knowledge  of  the  channel  feedback  history  is  then  asynchronous.  The 
objective  in  this  case  is  to  prevent  new  arrivals  from  interfering  with  some  collision  resolution  in 
progress.  This  is  possible,  if  each  user  can  decide  whether  a  collision  resolution  is  in  progress  or 
not,  within  a  finite  number  of  slots  from  the  time  he  generates  a  new  packet.  We  observe  that  a 
user  who  has  a  new  packet  and  observes  a  C  slot  decides  to  wait,  since  he  can  then  deduce  that 
there  is  some  collision  resolution  in  progress.  Also,  as  we  explained  in  Section  II,  the 
occurrence  of  two  consecutive  NC  slots  corresponds  to  either  two  consecutive  unit-length  CRIs 
or  to  the  end  of  a  CRI  which  started  with  a  collision  slot;  thus,  upon  the  observation  of  such  an 
event,  a  user  can  decide  that  there  is  no  collision  resolution  in  progress.  In  view  of  the  above,  we 


conclude  that  under  limited  feedback  sensing  the  algorithm  can  be  modified  to  operate  as  fol¬ 
lows: 

The  window  size  is  the  same  as  in  the  full  feedback  sensing  case.  However,  the  window  slides 
through  the  unexamined  interval  from  present  to  past  and  its  edge  is  maintained  one  slot  before 
the  current  time,  (see  Figure  2  ).  Within  each  window,  the  operations  of  the  algorithm  are  the 
same  as  in  the  full  feedback  sensing  case. 

In  the  limited  sensing  environment,  and  for  very  light  input  traffic,  the  algorithm  induces 
mean  packet  delay  equal  to  2.5.  As  the  rate  of  the  input  traffic  increases,  the  mean  delays 
approach  those  induced  under  full  feedback  sensing.  The  throughput  of  the  algorithm  remains 
identical  to  that  under  full  feedback  sensing.  In  Figure  2,  we  plot  the  expected  delays  that  the 
algorithm  induces,  under  both  full  and  limited  feedback  sensing.  In  the  latter  case,  the  mean 
delays  were  computed  by  using  methodologies  similar  to  those  in  [4],  [9], 


Remarks  We  point  out  that  the  modification  of  Capetanakis’s  dynamic  algorithm,  for  operations 
under  limited  feedback  sensing,  is  still  an  open  and  complex  problem.  In  contrast,  the  same 
modification  is  simple  when  the  proposed  algorithm  is  adopted.  In  systems  where  the  Poisson 
user  model  is  valid,  the  Part-and-Try  algorithm  with  binary  feedback  in  [8],  can  be  modified  to 
operate  under  limited  feedback  sensing,  [9],  [10].  The  throughput  of  the  latter  algorithm  is  then 
0.45.  But  when  the  Poisson  user  model  is  not  valid,  it  leads  to  deadlocks.  In  contrast,  the  pro¬ 
posed  algorithm  does  not  lead  to  deadlocks  and  its  performance  can  be  then  analytically 
evaluated. 


V.  The  Output  Traffic  Interdeparture  Distribution 

In  this  section,  we  concentrate  on  the  computation  of  the  output  traffic  interdeparture  distri¬ 
bution.  In  particular,  we  find  analytically  the  steady-state  distribution  of  the  distance  between 
two  consecutive,  successful  transmissions,  when  the  algorithm  in  this  paper  is  deployed.  We 
point  out  that  the  algorithm  generates  an  output  traffic  process  with  memory.  Our  computations 
correspond  thus  to  the  first  order  distribution  from  this  process.  This  first  order  distribution, 
together  with  a  memoryless  assumption,  can  be  used  as  an  approximation  of  the  actual  output 
traffic  process,  when  studies  of  systems  which  deploy  the  algorithm  and  interact  with  each  other 
are  undertaken.  Such  interactions  may  correspond,  for  example,  to:  (1)  Servicing  the  output 
traffic  from  several  systems  that  deploy  the  algorithm,  by  a  single  server  queue.  (2)  Transmit¬ 
ting  the  output  traffic  from  a  system  that  deploys  the  algorithm,  through  the  transmission  chan¬ 
nel  of  another  random  access  system,  (multi-hop  problem). 

Our  methodology  utilizes  the  regenerative  character  of  the  output  traffic  process  that  the 
algorithm  generates,  and  its  steps  are  as  those  in  [4].  We  define  the  sequence  {Prfm  as  follows: 
Each  P i  is  a  collision  resolution  point,  (CRP),  which  follows  a  slot  containing  a  successful 
transmission  and  at  which  the  lag  equals  one.  P\  is  the  first  after  zero  such  CRP,  and  for  every 
i>l,  Pi+ 1  is  the  first  after  P,  such  CRP.  Let  S(j  iSl  denote  the  number  of  successful  transmissions 
in  (0,  P^,  and  let  dn  denote  the  distance  between  the  (n-l)-th  and  the  n-th  successful  transmis¬ 
sion.  Then,  S„  1  is  a  renewal  process,  and  the  process  dn%  nSl  is  regenerative  with  respect  to  it. 


Let  us  define,  C/=Si+l  -Sit  »'£  1.  Then  C,  denotes  the  number  of  successful  transmissions  in  the 
interval  (/*,,  where  this  interval  will  be  called  the  i-th  cycle.  Let  us  define, 

f  1  ,  if  dH  =  s 

Ai(s)  -j  0  ,  otherwise  ^ 

H  =  E{Pi+l-PJ  (6) 

From  the  regenerative  theorem  [4],  we  then  conclude  that  if  C  =  E{C i }<  °°,  then, 

lim  AT1  £  In(s)  =  iim  AT1  Ef  £  /„(*);  =  C'1  £  /B(s);  (7) 

N~*~  n=l  <i=l  <i-l 

;  where  if  the  intensity  of  the  input  Poisson  traffic  is  X,  then, 

C  =  lH  (8) 

In  addition,  since  P(CX  =  1)  >  0,  the  distribution  of  C\  is  aperiodic  and  there  exists  a  random 
variable  d„,  such  that  the  sequence  dH,  n=l,2,...  converges  in  distribution  to  d».  Then, 
represents  the  steady  state  interdeparture  distance  induced  by  the  algorithm,  and  its  distribution 
satisfies  the  equality, 

.  c- 

P(d„=s)  =  C-'  EfZUs)}  (9) 

«-i 

c, 

The  finiteness  and  the  computation  of  the  quantities  C  and  £f£  /„(s)j  in  (9)  are  related  to  the 

»=i 

existence  and  computation  of  appropriate  solutions  to  infinite-dimensionality  linear  systems. 
Those  systems  and  their  solutions  are  included  in  the  Appendix.  In  Table  5,  we  include  the  com¬ 
puted  upper  and  lower  bounds,  respectively  denoted  P“  and  P‘s,  on  the  probability  P(d„  =s),  for 
input  traffic  Poisson  intensities  X  =  0.1  and  X  =  0.4.  In  Figure  3,  we  plot  the  lower  bounds  against 
s,  for  various  input  traffic  Poisson  intensities  X. 


s 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


X=0.1 

X=0.4 

Pls 

PI 

pi 
r  s 

Ps 

EBS9 

0.4702 

0.2000 

0.2048 

0.10806 

0.0641 

0.0696 

0.0787 

0.0537 

0.0603 

0.0603 

0.0502 

0.0591 

0.0280 

0.0397 

0.0420 

0.0503 

0.0196 

0.0264 

0.0364 

0.0452 

0.0099 

0.0173 

0.0332 

0.0431 

0.0057 

0.0094 

0.0265 

0.0393 

0.0013 

0.0062 

Table  5 


Upper  and  lower  bounds  on  the  interdeparture  distribution  for  Poisson  rates 

X  =0.1  and  X  =0.4 


From  the  results  in  Table  5  and  Figure  3,  we  draw  the  following  conclusions: 

(1)  For  low  rates  X  of  the  Poisson  input  traffic,  (X  <  0.1),  the  interdeparture  distribution  is  close 

to  the  Bernoulli  distribution  whose  parameter  is  p  =  Xe~x.  In  particular,  denoting  P,  ^  P(d„  =  s), 

we  then  have  P,~p  ,  for  s  >  2.  The  probability  P i ,  however,  is  then  significantly  larger 

than  the  Bernoulli  parameter  p.  The  intuitive  explanation  of  the  latter  phenomenon  goes  as  fol¬ 
lows:  For  small  rates  X,  single  arrivals  in  two  consecutive  slots  occur  with  probability 
p2=(Xe~x)2  =>  X2,  while  the  probability  of  a  collision  slot  is  then  approximately  equal  to 
2~l  X2  e~x  ~  2~l  X2.  Thus,  for  small  rates  X,  single  arrivals  in  two  consecutive  slots  contribute 
two  thirds-of  P\,  while  the  remaining  one  third  is  due  to  consecutive  departures  at  the  end  of  a 
collision  resolution  interval. 

(2)  As  the  rate  X  of  the  Poisson  input  traffic  increases,  the  interdeparture  distribution  induced  by 

the  algorithm  deviates  further  from  the  Bernoulli  distribution.  In  fact,  as  X  increases,  the  mass  of 

the  interdeparture  distribution  accumulates  at  relatively  small  s  values.  For  example,  for  X  =  0.4, 

10 

we  have  Px  =  0.471  and  £  Ps  ~  1. 

j-i 


Remarks  Our  results  show  that  it  is  generally  wrong  to  conjecture  Bernoulli  interdeparture  distri¬ 
bution.  In  fact,  the  obtained  distribution  is  far  from  Bernoulli.  Even  for  small  input  Poisson  rates 
X,  the  probability  P\  does  not  equal  to  Xe~x.  We  point  out  that  the  proposed  method  could  be 
applied  to  other  algorithms  as  well,  including  the  Capetanakis’s  algorithm.  However  the 
development,  of  the  corresponding  recursions  and  infinite  dimensionality  linear  systems  of  equa¬ 
tions,  is  then  very  complicated.  The  simple  operations  of  the  proposed  algorithm  present  an 
advantage,  which  does  not  characterize  other  existing  algorithms. 


10 


VI.  Conclusions 


We  presented  a  simple  window  random  access  algorithm  for  systems  with  binary,  collision 
versus  noncollision,  feedback.  We  analyzed  the  algorithm  in  the  presence  of  the  Poisson  user 
model,  and  under  both  full  feedback  sensing  and  limited  feedback  sensing.  In  addition  to  the 
throughput  and  the  delay  analyses,  we  studied  the  effect  of  feedback  errors  on  the  throughput  of 
the  algorithm.  We  also  studied  the  output  traffic  interdeparture  distribution,  (The  last  two  stu¬ 
dies  assume  continuous  feedback  sensing).  As  compared  to  the  Capetanakis’s  dynamic  algo¬ 
rithm,  the  proposed  algorithm  is  superior  in  terms  of  delays  and  insensitivity  to  feedback  errors. 
In  contrast  to  the  former,  the  algorithm  can  also  be  easily  adapted  for  implementation  under  lim¬ 
ited  feedback  sensing,  it  allows  for  analytical  studying  of  the  output  traffic  interdeparture  distri¬ 
bution,  and  can  be  easily  analyzed  when  strici  delay  limitations  are  imposed,  [3]. 


VII  Appendix 

VII.  1.  Stability  Analysis 

We  will  present  the  stability  analysis  in  the  case  feedback  errors  may  occur.  By  setting  e  =  5  =  0 
in  our  results,  we  get  the  corresponding  quantities  in  the  error  free  case. 

The  stability  region  of  the  algorithm  is  provided  by  inequality  (1),  where  the  expected 
value  in  it,  is  given  by  expression  (2).  We  start  with  the  computation  of  the  expected  values 
Lk  ;  k  >  0  in  (2a). 


Computation  of  Lt 
We  define: 


Gn  k-n'  The  expected  number  of  slots  needed  by  the  algorithm,  for  the  successful  transmis¬ 
sion  of  k  packets,  given  that  n  of  those  packets  have  counter  values  equal  to  one,  and  k-n  of 
the  packets  have  counter  values  equal  to  two. 


Notice  that  Lk=G^0,  for  k>2,  while  Lk  * G^ Q,  for  k  <2.  We  first  show  how  to  compute  LQ 
andLi. 


Computation  of  Ln: 

From  the  operation  of  the  algorithm  we  have 

1 ;  w.p.  (1— e) 


where 


From  (Al)  and  (A2)  we  find  that 


L°  1  1  +G0,o  ;  w.p.  e 


'1  +L0  ;  wp.  (1-e) 
1  +(70,o  •  WP-  e 


Ln  = 


(1-e)2 


(Al) 


(A2) 


(A3) 


Computation  of  L  i 

It  was  found  that  L\  satisfies  the  following 


Lx 


1 ;  w.p.  (1-5) 
I+Gio ;  WP-  0.58 
l-K/o.i  ;  w P-  0.58 


(A4) 


where  G  i,o  and  G0.i  satisfy  the  following 


1+L0  ;  w.p.  (1-5) 
Gjfo  =•  l+G i_o  ;  w.p.  0.55 
1+C/o.i  :  w.p.  0.55 

» 

'1+Lj  ;  w.p.  (1 — e) 
<^0-1  1+G0,i ;  w.p.  e 


From  (A4),  (A5),  and  (A6)  we  find  that 


1+  — - —  +  — § —  i+-M- 
2(l-e)  (2-5)  +  (i_e>* 


5 

2(l-e) 


Computation  of  U.  for  k>2. 

From  the  operation  of  the  algorithm  we  obtain: 


G„,k-n  =  1  +  1  J  r*  G,-*.,  :  n>2,k>n 

i=0  L  J 


where 


Go-*=a^r+L* 


G<  *  = 


l+L*(l-6)  +  0.55  -7T—+Lk+l 
(l-€) 

1-0.55 


It  can  be  shown  by  induction  that  Gn  k_n  has  the  following  form 

^/i,*-(i=A^G*,o+A?)<7ii_i>o+A5i5)  ;  2<n<k 
;  where  a£‘\  i=l  ,2,3  are  independent  of  k  and  can  be  computed  recursively  as  follows: 

Ac D_J±§_  4(1-5)  m  1 4-35-1 2e  +  4e5 

*  '  6-36  ’  *  "  6-35  *  *  3(1— e)  (2-5) 


^.[i-2-r'r1-  .  m3 

4®  =  f  1-2“"1-1  2" "«  +  y  [n]  4(2)1  _> 3 

122  \  (1-0.55)  itS  l  *J  ‘  J  ’ 


=  [l-2-"l-1  1+  — — - +  — - + — - +2~n  V  fn]  n> 3  (AIS) 

*  1  1  [  (1-0.56)  (1-e)  (2-6)(l-e)  "  U'l  '  J’  (A15) 

For  n=k,  expression  (All)  gives: 

A?'  a<?) 

Lk-Gk, 0-  *•*-'  +  (A16) 

Expression  (A  16)  together  with  the  recursions  in  (A13)-(A15)  provide  a  mean  for  the 


(A10) 


(All) 


(A  12) 


(A  13) 


(A  14) 


(A  15) 


(A  16) 


computation  of  Lk,  k22. 

Development  of  an  upper  bound  on  L, 


It  can  be  seen  by  induction  that: 


4”^ -2^  :**2 

(A  17) 

=  I  ;  *22 

i-4‘> 

(A  18) 

a13)s  (2_2e+5)  k  +  1 
*  (2— 5)(1— e)  (1— e) 

(A19) 

From  (A16)-(A19)  we  finally  find: 


,  <  /  ,  3(2-2e+5)  3(2-6) 

*  *"*  (l-e)(4-58)  *  (1— e)(4— 55) 


,  k>2 


(A20) 


from  which  we  can  show  that 


LkZ 


3(2-2e+8)  fc2  +  __3(6-2e-5) 


2(1 -€>(4-5 5)  2(l-e)(4-58) 


it  + 


L,- 


6(2— e) 


(l-e)(4— 55) 


^Lj  =  aA2+p*-Hy  (A21) 


For  e  =  5  =  0,  we  obtain  the  quadratic  upper  bound.  Due  to  the  upper  bound  on  Lk,  we  conclude 
then  that  the  following  condition  is  sufficient  for  stability: 

30 

£  Lk  p(k  I  A)  +  £  Li  p{k  I  A)  <  A  (A22) 


*-o 


*-31 


where 


p(k\A)-e 


-u  (XA)* 


A:! 


After  some  manipulations,  we  conclude  that  (A22)  is  equivalent  to: 


30 

/(AA)=£L*p(/klA)  +  a^ 

*. o 


„  30 

(XA)2  +  \A-  £  X2p(*IAn  + 
*-o 


+  P 


30 


XA-  £  *  p(*)A)| 
*=o 


30 

1-  £p(*!A) 
*» 0 


<  A 


(A23) 


Let  us  now  define: 


*  AXA 


(A24) 


Then,  from  (A23)-(A24)  we  conclude  that,  for  the  stability  of  the  algorithm,  it  is  sufficient  that 
the  input  rate  X  satisfies  the  following  inequality. 


(A25) 


\<sup 


X 

fix) 


The  following  condition  specifies  a  region  of  X  values  for  which  the  algorithm  is  unstable. 


30 

where  #(x)  =  £ 

k=0 


Lk  e~x  xk/k\ 


X  >  sup 
*20 


X 

six) 


(A26) 


The  maximization  of  expessions  in  (A25)-(A26)  has  been  done  numerically,  and  provides  the 
throughput,  as  well  as  the  optimal  window  size  A*.  In  all  the  cases  the  order  of  the  difference 
between  the  two  suprema  in  (A25)-(A26)  is  less  than  10-3.  The  optimal  window  size  is  found  as 
x *  (X*)"1,  where  x*  is  the  value  that  attains  the  suprema  in  (A25). 


Computation  of  L\.  for  the  Capetanakis’s  Dynamic  Algorithm 
Here,  the  quantitites  Lk,  k>2  can  be  computed  recursively  as  follows: 


L\  =  [1-2  2-*]-1 1  1+2  Lq  2“*  +  2 


*-i 


f*l 

*> 

2“* 

I'J 

* 

where 


.  •  1  l-2e  +  5 

Lo  ~  7TTT  L\  =■ 


(l-2e)  *  (1— 2e)(l— 6) 

Moreover  the  following  upper  bound  on  has  been  found. 


Lk< 


3(1— e)  (  2(5-e) 


(1— 2e)  (1— 2e)(l— 6) 


k-  1 


(A27) 


(A28) 


(A29) 


Using  the  bounds  in  (A29),  we  computed  upper  and  lower  bounds  on  the  throughput  for  the 
Capetanakis’s  dynamic  algorithm.  The  computed  upper  and  lower  bounds  were  identical  to  each 
other  up  to  the  fourth  decimal  point,  and  are  included  in  Table  3. 


VII.  Delay  Analysis 

The  following  definitions  will  be  used  in  the  sequel. 

/  :  Length  of  a  conflict  resolution  interval 

x  :  Window  size. 

E(X/x )  :  Expected  value  of  the  random  variable  X,  given  that  the  window  size  is  x. 

p  (x/x)  :  The  probability  that  the  conflict  resolution  interval  has  length  x  given  that  the 

window  size  is  x. 

gd  :  Number  of  slots  needed  to  reach  a  CRP  with  lag  1  given  that  the  current  lag  is 

equal  to  d,  dcF. 
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wd  :  Cumulative  delay  experienced  by  all  the  packets  that  were  successfully  transmit¬ 

ted  during  gd  slots. 

N  :  Number  of  packets  transmitted  in  the  CRI,  that  starts  at  time  t. 

z  :  Cumulative  delay  of  the  N  packets,  after  the  CRP  t. 

V  :  Cumulative  delay  of  the  N  packets,  until  the  instant  t2- 

Let  us  also  define, 


Gd  =  E(£d)  Wd  =  E  (wd) 

Note  that  by  definition, 

g. 

Wl=E(£Dn)  =  W 

nm  \ 

Also,  the  mean  session  length  G=E(Ti+l-Ti) ;  i  £  I,  equals  to  Gx.  If  G{  <  then  by  Wald’s 
identity  we  have  that, 

<2=XG, 

Therefore,  the  determination  of  W\  and  G\  will  permit  the  computation  of  the  mean  packet 
delay. 

The  operation  of  the  algorithm  yields  the  following  relations  for  the  gj,,  d  eF. 


1  £  d  £  A  gd  =  l  if  /  =  1 
lSdSA  gd  =  l+gt  if  /  >  1 

and  that 


A  <d  gd  =  l  +  gd-*+i 
Taking  expectations  in  (B1)-(B2)  yields: 

Gd=E(l ld)  +  £  G,p(s/d )  \i  \<d<A,  deF 

t*\,ssF 


(Bl) 

(Bla) 

(B2) 

(B3) 


Gd  =  £(/ 1  A)  +  Y,Gd-tk  +  sP  (* /A)  if  d>  A,  deF 

JCF 


(B4) 


Equations  (B3)  and  (B4)  comprise  a  denumerable  system  of  linear  equations.  Of  interest  to  us  is 
the  element  Gi  of  a  particular  solution  of  this  system.  We  now  proceed  in  the  development  of  an 
initial  upper  bound  on  the  solution  of  the  system  in  (B3)-(B4).  Following  the  methodology  in 
[4],  such  a  bound  will  be  the  sequence  Gdu  =  yud+Z,u,  if  yu,  can  be  determined  so  that  the  fol¬ 
lowing  inequalities  are  satisfied 


G°d'UZE(l\d)+  £  G°,up{sld)  =  Gdu 

s*\,stF 


\<>d<A,  deF 


G%u>E(tlA)+  ZG2_a  +  jP(s/A)  =  G}'U  if  d>A,  deF 


seF 


(B5) 


(B6) 


Substituting  Gd  in  the  right  hand  side  of  inequalities  (B5)  and  (B6).  it  can  be  seen  that  if  deF, 


Gi*  =  G%u  +  y*(£(//d) -d-  (1 +U)e~u)  -  ^(l+Xd)e~u  if  \<d<A 
Giu  =  Gd'U  +  E (l  I  A)  -  Y«  (A  -  £  (//  A))  if  d> A 


(B7) 


(B8) 


From  (B8)  we  conclude  that  if  £(// A)  <  A,  the  condition  for  stability  of  the  system,  inequalities 
(B6)  are  satisfied  if. 


£(/ia) 

1  A  -  £  (//A) 


(B9) 


With  this  value  of  yu,  it  can  be  seen  that  inequalities  (B5)  are  satisfied  if 

C«  =  ma xf-yu,  sujp(Q(d))} 


where 


ec d)  = 


£(//  d)  +  yu(E(l\d)-d-(\+\d)e-u] 
(\+ld)e~u 


(BIO) 

(Bll) 


From  the  above  discussion  we  conclude  that  the  solution  to  system  (B3)-(B4)  satisfies  the  ine¬ 
qualities 

GdZyud  +  Zu,  dtF  (B12) 


The  uniqueness  of  the  solution  is  guaranteed  by  the.  same  techniques  as  in  [4],  If  we  use  a  simi¬ 
lar  method  for  the  development  of  a  lower  bound,  we  find  that 

Y/4  +  Ci  =  G,,4eF  (B13) 

where 

Jt  =Y«,  and  bm  (B14) 

From  the  operation  of  the  algorithm  we  also  have  that  Wd ;  dz F  satisfy  the  following  system 
of  linear  equations. 


Wd  =  E(z  \d)  +  E(y\d)+  £  Wsp(s\d),  l<d:£A,  dz F 

seF,j  *  I 

tVd  =  £(2lA)  +  £(vlA)+2:iV(<_i  +  ,p(slA),  A  <d,  dz F 

jeF 

Following  the  methodology  in  [4]  we  can  show  that 

W°dJ  =  p/d2  +  V/d  +  %,<>Wd<,  pud2  +  vud  +  §„  =  W$H 

where, 

XA 


Hu  =  Uf  = 


2(A  -  £  (I  t  A)) 


(B15) 

(B16) 

(B17) 

(B18) 
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m  JL*  AJI  AA.' 


(B 19) 


v«  =  V/  = 


£  (z  I  A)  +  £  (y  I  A)  -  XA2  +  |iu£  ((A  -  /)2 1  Aj 
A  -  £(/ 1  A) 

=  sup  (<K<0),  £/  =  in/-  (<>(<*)) 

1  £  ok  A  1 £  a  £  A 


(B20) 


where. 


E (2  I 4)  +  £ (V I 4)  +  (l 2 1 d)  -  d 1  -  ( 1  +  \d)e-“)  -  v„(4  -  E  (/ 1  4)  -  ( 1  +  U)e "“)  „ 

<t>(a)  = - ~ — TTT-u -  (B2 1 ) 

(1  +  A4)«  ** 

Substituting  in  the  right  hand  side  of  (B  15)-(B16),  we  obtain  W°  .  We  find, 
tVl,„  =  E(z  1 1)  +  E(y  1 1)  +  p*(£(/2 1  lHt  +  X)«"X)  +  vm(£(;  1 1)  -  (1  +  X)e'x)  +  ^(1  -  (1  +  \yx)  (B22) 

w\ J  =  W\'U  -  (%„  -  %,)  (1  -  (1  +  X)e'x)  (B23) 


Also, 

G\m  =  £ (/ 1 1)  +  Ym(£ (/ 1 1)  -  (1  +  X)e_x)  +  c«(i  -  a  +  X)<?_x) 

ci./ =<?!.„ -(C«-&)0- a +*)*'x) 

From  the  regenerative  theorem  [4],  we  have  that 


,  W,  w  j  u 

D;  =  — -i—  S  D  £  — —  =  r»“ 


W\ 


to 


XG* 


=  D“ 


(B24) 

(B25) 

(B26) 


l,/ 


In  this  Appendix  we  also  show  that  the  conditional  expectations  of  the  form  E(X\d)  can  be 
computed  with  high  accuracy.  Let  us  define 

£ (X  I  d,k)  :  The  conditional  expectation  of  the  random  variable  X,  given  that  the  arrival  interval 
contains  k  packets,  and  has  length  d.  Then, 

E(X\d)=  JjE{X\d,k)e~u^jL  (B27) 

The  quantities  E(X\d,k)  depend  only  on  k.  In  Section  V.I  we  show  that  the  quantities 
£*  =  £(/ 1 d,k)  can  be  computed  recursively: 

£0  =  L  j  =  1  (B28) 


4  = 


Aj?> 

1  -ApLk~' 


+ 


_AP_ 

i-Ay> 


;kZ2 


where  a£),  «=1,2,3  can  be  computed  recursively  as  follows: 


(B29) 


n>  3 


(B30) 
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A^=[l~2- 


(B32) 


"]j  l+T*(n+iy*T*  2  A<3>Rk  n> 3 

(i=2)  L  JJ 


a£»=1/3.  A?>  =  2/3,  A?>  =  7/3 

The  quantities  Zk=E(z\  d,k),  can  be  also  computed  recursively: 

Zo  =0,  Z\  —  1 

45>  _  aV 


z>  = 


-Zk- 1  + 


:  kt 2 


‘  l-A^ 

and  the  quantities  a£\  i  =4,5,6  can  be  computed  recursively  as  follows: 
a<,4)  ;  n  >  3  satisfies  the  recursion  in  (B30) 

Ai5) ;  n  >  3  satisfies  the  recursion  in  (B31) 


A?> 


=  [l-2_,,r,n  +  n2-"  +  rt22"’,+2-',£^  A^}  ,  n>  3 


(B33) 

(B34) 

(B35) 


(B36) 


44>  =  f  A?>  =  Af  =  \ 


The  quantities  Yk  =E(l2  \d,k)  can  be  computed  as  follows: 


Yf 


Y0  =  0.  Yi  =  1 

AP  „  A?) 


-ft-i  + 


;  k>2 


1-Aj?  *-*  '  !_aP 

and  the  quantities  A^  ,  <  =  7,  8, 9  can  be  computed  recursively  as  follows: 

A^  ;  n  >  3  satisfies  the  recursion  (B31) 
aJ,8)  ;  n  >  3  satisfies  the  recursion  (B32) 

A?  =  [1  -  2-T72A.-1  +  2-(l  +  2L„)  +  «2-',(l  +  +  2-'I'£  A< \9)}  n  >  3 


(B37) 

(B38) 

(B39) 


(B40) 


A?  =  A?>  =  |,  Af  =  16  (B41) 

From  the  above  formulas,  we  see  that  a  finite  number,  M,  of  terms  from  the  infinite  series 
in  (B27),  can  be  computed.  Also,  for  large  k  values,  and  based  on  the  recursive  expressions,  sim¬ 
ple  upper  and  lower  bounds  on  E(X  I  d,k)  can  be  developed.  Those  bounds  can  be  used  to  tightly 
bound  the  sum 

£  E(X\d,k)e~u^- 

k-M+ 1  *• 


Remark  It  can  be  also  proved  that  E (y  I  d)  =  ■ 


H1* 


V.ni  Interdeparture  Distribution  Analysis 
We  first  provide  some  definitions. 

lk_m:  Given  k  packets  with  counter  values  equal  to  1  and  m  packets  with  counter  values 

equal  to  2,  the  number  of  slots  needed  by  the  algorithm  until  the  first  successful 
transmission,  (and  including  it),  after  the  k-multiplicity  collision  has  been 
observed. 

Given  a  collision  resolution  interval  which  starts  with  a  k-multiplicity  collision, 
the  number  of  length  s  interdeparture  intervals  within  it.  The  length  from  the  ini¬ 
tial  collision  to  the  first  successful  transmission  is  included  in  the  counting. 

Starting  with  a  CRP  at  which  the  lag  equals  d,  d>l,  and  which  follows  a  success¬ 
ful  transmission,  the  number  of  slots  needed  by  the  algorithm  to  reach  the  first 
lag-one  CRP  which  follows  a  slot  containing  a  successful  transmission. 

Starting  with  a  CRP  at  which  the  lag  equals  d,  d>  1,  and  which  follows  a  success¬ 
ful  transmission,  the  number  of  length  s  interdeparture  intervals  until  the  first 
lag-one  CRP  which  follows  a  slot  containing  a  successful  transmission.  The  dis¬ 
tance  from  the  initial  CRP  to  the  first  successful  transmission  is  included  in  the 
counting. 

Given  an  arrival  interval  of  length  d,  the  probability  that  there  are  k  arrivals  in  it, 
that  /*,o=S.  and  that  it  takes  /  slots  for  its  resolution,  including  the  initial  collision 
slot. 

/’*(/):  Given  a  k-multiplicity  initial  collision,  the  probability  that  it  takes 

/  slots  for  its  resolution,  including  the  initial  collision  slot. 

The  above  definitions  are  needed  for  the  derivation  of  recursions  that  are  pertinent  to  the 
infinite-dimensionality  systems  associated  with  the  quantities  in  (9).  We  first  note  that: 

H  =E{h\}  ,  C  =  IH 

c, 

n= 1 


nky 

hd\ 

mdy 

P(*./.  8  Id): 


Auxiliary  Recursions 

The  operations  of  the  algorithm  induce  the  following  recursions: 
(I)  l].m-0,Ym  ,  P(lkm=0)=0  ,Y k>2,  Vm 

f  1  ,  if  m= 1 

to.m  =  1H*,0  .  P(h,m  =  1)  =j  0  ,  \{m*\ 
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'  A 

km  =  1  +  km+k-i  :  probability  K.  2~k,  k>2 

0  - 


1  ,  if£=l  and  5=0 

P(km=s)=\  P(lm,o  =  s~\),  if A:=0,  5>1 


2~*  £  { P  (k  m+lt-i  =  5-1),  if  *>2,  5  —  1 

i=0  k  J 


f  1  ,  if5=l 
n,*'=  1  0,  if  5*1 


•  w/r/j 

£^2  ;  «*  ,=  , 

*  1  +nk-i,s  ,  vv/i 


.  with  probability  P(ko*s-l) 

1  +nk-i  s  ,  wu/i  probability  P(k,o  =5-1) 


£  P(l,,o=s-l)  ,  if5>l 
A  *=2 

Nk.s  t  E  K  i  * 

l  +  £/5tt.O=0)=l,  if 5  =  1 

1=2 


fl  .  if /  =  1 

*=°,  1  ;  />*(/)=  0  f  otherwise  ,  P2(l)  =  Pill o  =  1-2),  for  />3 


it >2  ]  /-*-i 

/>*+2f:/>*(/)=  £  ^M>=5)/V,  (/-5-1) 


Given  Poisson  intensity  X, 


P(U,p\d)=  e~*  1jf~P(ko  =  P)Pk-i  (/-P-1) 


Recursions  for  h,, 


J  <  I  i  1  ,  (  i  I  i  "  » J.t  .  ■ 


«  "**  V  V. 


r,  i’  i*  *%/*#*»>  1 

*  *  St  4:)1  ?  1,-  * 


=*1 


V, 


*  (f»  .  (  , 


Given  Poisson  intensity  X,  and  from  the  operations  of  the  algorithm,  we  easily  conclude: 


1  ,  with  probability  hie  ^ 
d<A  \  hd  -  <  1+A  i  ,  with  probability  e_XA 

l+ht ,  with  probability  Pk(l),  l> 2 

v  k=2  k\ 


d>A  ;  hd  =  l+hd_A+l  ,  with  probability  u  Pk(l) 

*to  *! 


and  thus. 


Hd^E(hd}=  e~u  +  E{l\d}+e-uHl  +  'Z  Ze~U  ^TT  P^H‘ ; 

*>2  />2  *• 


Hd=  E{l\A}+£  £e 
*20 />i 


-U  (XA)* 


Ek(l)Hd-\+i  d  >  A 


where. 


I  £  e-u®4tpkW 

lrX\  r>\  k  • 


Recursions  for  md , 

For  [J  denoting  integer  part  ,  for  w.p.  meaning  with  probability,  and  for  Poisson  intensiy  X, 
we  conclude: 

r 

«i,j :  w.p.  e~^Xd 

nk,s  +  mi,  ;  w.p.  e~u^~-Pk{l)  ;  k> 2 

I 

dZA\md',=  nk-hs+mu  \wp.P (01  d)  n>0  Pn(0ll)P(k,l,p  I  1) ;  k>2 

n+p+2*s 

,  i*  I 

1  +«*_!,,  ;  w.p.  e-Md+])—  n> 0  <rX',/>(4,0=p)P*_l(/-p-l);X>2 

*  ■  n +p+2^s 

.1  ;  w.p.  Xe-hd+s-i)  .  5>2 

For  d  >  A: 


=  nk'S+md-M.s  ;  w.p.e"^  — r~  Pk(l)  ;  *>1 

! 


~  nk-),s  +  WP-  Z  e 


-XA(n -t-1)  1^.  P(lk  Q=p)  Pk_{(l-p-\)'  lc>lt  if  L-— -J  > 


'l+nKi  +  ^U-iW*  :wP-  I  e-^(^1)iMli/>(/fc0  =  5-rt-l)P*_I  (/-J+/I),  k>\, 

d~A  *  • 


«*-u  +m,s  ;  w.p.  e 


4^1  w-i^fj(i-D 


:*rL  A=fJ  21 


Z  P(/*.o  =  P)/>*-i(/-p-D;^>2 


l~ji  +p+i« 


=  1  +  «*-!, x  +  mu  ;  w.p.  e 


a-i  /-m44j  ;/t52 

A-l 


«*-u  +  m.s ; 


\  ,  i  d  —A  i  . . 

-*4H-r-rd  +l) 

A_1  n>0 


,-x*  X* 


'*+2+p-+i~f'J  « 


^  7r/>(/^o  =  P)^-1(/-p-i);^ 


_X-  +1  X*  f  ii—A 

\ +nk.Us+mu  iw.p.e  Z  «“*"  77  ^  4,0  =J-«-2-[  — H  . 

fi^O  *•  l.  A_1  . 


A-i  /+"-^+i:~7-J+1  :**2 


=  1  ;  w.p.\e-X(d+s-"  ;  if  j— 2-1  ^-=-J  >0 


■  u«  iiTwiMwiru in«n 


# 


L 


f 


f» 


4* 


* 


C 


=  1  ;  w.p.  X 


d< 


d-b 


A-l 


Let  us  define. 


A\\  ,  x>0 

"»=  0.*<0 


PM  i  I  e-*  77  P»(/) 

*21  *  • 

*8.4  S  **.,  ^  T7  •  *8.1  =  l-«- 

*21  *  1 


A  r-,  ,-s  5* 


I  *“*^7  ^4.0  =  5) 


*21 


Then,  using  the  above  defined  quantities,  and  the  recursions  on/nd 


For  d  <  A: 


A  e~u 

Md.s  t  E{md.s}=  N  u,  +  X  A*u  l^uCD  +  fr 

/23  l-« 


Px  (01+ 


l-e 


-x 


Nx.s-P\.s- 1 


f-2 


+  l/(j-2)e~*4|  g-xo-2)  g 


m=0 


+  -  X 


l-e' 


For  d>  A  and  p  ^  |_  ^~|  ,p=0,l,.„: 

Mj.i  =  *U,J  +  X  Md-M,s  P xa (0 
/21 


+  g  J  *X(d+p-pA),j  P\(d+p-p&),s-\  4"  £  P X(d  +p -pA)(0 

l  '23 


g  -XW+p) 
l-e“x 


*x.,  -Pk,-i  +  XMijPiW 


1 23 


we  easily  find: 


T  X.m  + 


4 
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^XA,j-^XA.j-lj  +  Z  ' 
L  J  i  <*<*> 

?-X(j-1) 

s-2 

m=5-l— min(p,  5—1) 


-  U(s- 2)\e~Md+p) 


( d+p-pA)e + 


l-e 


-x 


+  Uis-l^e^  Pu<t+p-ps).s-i- 


s  -2-p 

+  U(s-2-p)  e~Md+s-2)  £  e  p\j 

m  =0 


Bounds 


For  the  numbers  N^s,  we  used  the  following  bounds: 

0<Nk%s<k-\  ;  Ys 

Regarding  the  numbers  Hd,  we  used  the  methodology  in  [4],  and  proved  that. 


where. 


c lid  +  P/  <  Hd  <  0^  d  +  ,  d  >  1 


Ct^  =  au  =  (A-£:/r/IA;^,  E(l  I  a; 


P'  =  l<d<A2(d)  ’  P“  =max 
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Bounds  on  the  numbers  MdiS  can  be  developed  similarly  with  those  for  the  numbers  Hd.  The 
former  are  significantly  more  complicated,  however.  Instead,  we  used  the  following  simpler  and 
intuitively  clear  bounds,  where  Hd  denotes  the  upper  bound  on  the  quantity  Hd : 


0  <Mdj  3 


We  used  the  above  bounds,  for  d  >  30. 
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Figure  1 

Window  Selection  in  the  Limited  Sensing  Environment 


Lower  Bounds  of  the  Interdeparture  Distribution 
for  Various  Poisson  Rates 
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